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We study multilayer triangular lattice Ising antiferromagnets with interlayer interactions that are 
weak and frustrated in an abe stacking. By analysing a coupled height model description of these 
systems, we show that they exhibit a classical spin liquid regime at low temperature, in which both 
intralayer and interlayer correlations are strong but there is no long range order. Diffuse scattering 
in this regime is concentrated on a helix in reciprocal space, as observed for charge-ordering in the 
materials LuFe 2 04 and YbFe 2 04 . 
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Some simple models of highly frustrated magnets de¬ 
velop strong correlations at low temperature without 
long-range order. In this regime they are known as co¬ 
operative paramagnets or classical spin liquids [T]. The 
Ising model with nearest neighbour antiferromagnetic in¬ 
teractions on the triangular lattice was one of the earli¬ 
est examples to be studied in detail [2], while the Ising 
antiferromagnet on the pyrochlore lattice is an example 
of high current interest as a model for spin-ice materi¬ 
als [3]. The description of the cooperative paramagnetic 
state presents a theoretical challenge, and may involve 
emergent degrees of freedom and fractionalized excita¬ 
tions: a height [I] or gauge field [5] in these cases, with 
vortex or monopole excitations. 

In both these models the cooperative paramagnetic 
regimes are continuously connected to disordered ground 
states with macroscopic entropy. By contrast, systems 
that have ordered ground states typically do not dis¬ 
play strong correlations without long-range order, ex¬ 
cept near a critical point. Here we study a remarkable 
exception: the three-dimensional Ising antiferromagnet 
built from a6c-stacked triangular lattices, with interlayer 
coupling J_\_ much weaker than the in-plane coupling 
J. For = J this is the face-centered cubic lattice 
model, which orders discontinuously |5] at a temperature 
Tc ~ 1.74J, while for Jj^ = 0 it reduces to uncoupled 
layers, which remain disordered to T = 0. For J± J 
we show that there is a temperature window in which the 
model has strong correlations, both between and within 
layers, and a large but finite correlation length. We for¬ 
mulate a theory for this regime in terms of coupled height 
fields. It is striking for its correlations (helical in re¬ 
ciprocal space), fluctuations (quartic in wavevector) and 
mechanism for suppression of long-range order (bound 
vortex-antivortex pairs). It also has interesting parallels 
with theories for smectic liquid crystals |7] and for frus¬ 
trated quantum magnets [5]. 

To put this layered antiferromagnet in context, note 
that the ordering pattern in magnets is determined at 
mean-field level by the location of minima in the eigen¬ 
values of the exchange interaction matrix J(q) as a func- 



FIG. 1: (a) The multilayer triangular lattice with abc stack¬ 
ing. (b) The mapping (modulo 6) from ground state spin 
conhgurations on sublattices A, B and C to integer-valued 
height variables h (adapted from [2]). 

tion of wave vector q. The suppression of ordering in 
many frustrated models is a consequence highly degener¬ 
ate minima. As examples, nearest neigbour (nn) interac¬ 
tions on the pyrochlore lattice lead to a degenerate mini¬ 
mum band spanning the full, three-dimensional Brillouin 
zone [5], while competing first and second neighbour in¬ 
teractions on the diamond lattice give rise to minima on 
surfaces [S]. On the a&c-stacked triangular lattice, nn in¬ 
teractions have long been noted for generating minima on 
helical lines in reciprocal space [10] . For the Ising model 
on this lattice we show in the following that reciprocal 
space correlations are concentrated close to such a helix 
over an extended temperature range. 

We start from the model illustrated in Fig. , with 
the Hamiltonian 

H = J J- Tl 0'i,zO’j^z+l ■ (1) 

Here = ±1, the integer z labels planes, nn pairs of 
sites in the same plane are denoted by (ij), and those 
from different planes by {ij}- Of the possible further- 
neighbour interactions that we have omitted, the most 
important is the unfrustrated coupling 

Hs = Js cri,zO'i,z+3 

i,z 


( 2 ) 
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between a spin in layer z and the spin directly above it 
in layer z + 3, which lifts the degeneracy of the helical 
minima in J(q). 

To discuss the low-temperature behaviour of this 
model we exploit the fact that ground states of a sin¬ 
gle triangular layer can be described by a height model 
mm- Excitations out of the ground state are repre¬ 
sented by screw dislocations or ‘vortices’ in the height 
field. Vortex-antivortex pairs, present at any finite T, 
are unbound in an isolated layer and the vortex separa¬ 
tion f sets the correlation length, which is large at low T 
because they are dilute. 

Interlayer coupling greatly reduces the degeneracy. 
Ground states [12] for Jj_ ^ 0 are described by a height 
field with a gradient of maximum magnitude but arbi¬ 
trary direction. Of these, thermal fluctuations favour six 
states m in which each triangular layer contains alter¬ 
nating stripes of up and down spins. For low T the sys¬ 
tem adopts one of the three-dimensionally (3D) ordered, 
symmetry-breaking phases of this type. 

At intermediate temperatures, thermal fluctuations 
compete with the interlayer coupling. For the regime 
of most interest, J_l <C T ^ J, it is necessary to take 
two aspects of the physics into account. First, per¬ 
turbing around a system of uncoupled layers, we use a 
renormalization-group (RG) calculation [T3| to treat the 
couplings K± = J±/T and = J 3 /T, which are RG- 
relevant. Second, we make a detailed analysis of the 
influence of defects that appear in the coupled height 
model at finite temperature. The outcome depends on 
three quantities: the vortex separation the scale £± 
at which the renormalised coupling K± ^ 0{1) in the 
vortex-free system; and the strength of the renormalised 
K 3 . Interlayer coupling without vortices leads to 3D or¬ 
der at the scale £±. In addition, it generates a potential 
that confines vortex-antivortex pairs. In many settings, 
unbound vortices destroy long-range order but bound 
vortex-antivortex pairs do not. Remarkably, we find in 
this model that 3D order is disrupted even by bound 
pairs if K 3 is smaller than a critical value K^. Hence 
three distinct regimes of behaviour emerge from this anal¬ 
ysis, with: (a) weakly-coupled paramagnetic layers, when 


^ (b) 3D order, when .^_l <C ^ and K 3 > iF|; and 

(c) strong interlayer correlations but no long-range order, 
when £_l <C ^ and K 3 < K^. 

We now set out these calculations in more detail. For 
a single triangular layer, the mapping mm between 
ground-state spin configurations and single-valued height 
variables, defined at the centres of triangles, is illus¬ 
trated in Fig. Bb)- The inverse mapping has the form 
CTa = fih+Sa), where f{h+6) = f{h) with f{h) = -fl for 
h = —1,0,1 and f{h) = —1 for h = 2,3, 4. Here the sub¬ 
lattice label a = A, B or C and = 0, sb = 2, sq = —2. 
The dominant contribution to long-distance correlations 
involves the lowest Fourier component of f{h), so that 

TT 

aa ^ COS-{h + Sa). (3) 

Excitations out of the ground state consist of triangles 
in which all three spins have the same orientation, and 
the height field changes by ±6 around a closed path sur¬ 
rounding one such excitation. 

The probability of a coarse-grained height configura¬ 
tion h{r) for an isolated layer is proportional to 
with an effective Hamiltonian mm 

"^ 20 = J |y |V/i(r)|^ - gcos27r/i(r)| . (4) 

The first term in Eq. Q describes the entropic cost 
of height gradients; the second term encodes the fact 
that microscopically the height field is integer-valued. 
The value of the stiffness K is fixed by comparison 
with the known asymptotic behaviour (CTa(r)(T^(0)) ^ 
r“^/^cos|(sa — sp) of ground-state correlations in the 
Ising model US): for K = ^ this form is recovered and 
g is RG-irrelevant. Note that, because of the sublattice- 
dependent phase Sa in Eq. ([^, small wavevector fluc¬ 
tuations of h{r) represent Ising spin fluctuations with 
wavectors near the corners of the triangular-lattice Bril- 
louin zone. 

To describe a multilayer system we introduce a height 
field hz{v) in each layer, with r = {x,y). Using ([^, the 
inter-layer coupling written in terms of heights is fl6j 


/ ^2 

-y (dj;h^{r) cos ^(/i 2 +i(r) - ^^(r)) - sin ^(/i^+i(r) - h^{r)) + ...^ (5) 




where £ is the short-distance cut-off and the ellipsis indicates omitted terms that are RG-irrelevant at weak interlayer 
coupling. We introduce the notation | — ^^(r)) and define a reduced interlayer coupling ft:j_ = 

Ki_l£K. The leading contributions to the multilayer height Hamiltonian can then be combined as 

'Hzt! = j d^xiy{da:h,,(y) - K^cos5^P{v)^ -f -f Kj_sin(5(^)(r)^ -k 


( 6 ) 
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This model has a striking continuous symmetry under 
a joint real-space and height-space transformation. Let 
Re denote rotation in the xy plane by the angle 9. Then 
the transformation r' = i?e(r) and 

hzir) K{v) = hz{r') + Wz/n (7) 

leaves Hsu invariant for any 0. Ground states of Hsu are 
parameterised by constants 0 and a, having a form 

= k±{xcos9 — ysin9) + Wz/tt + a (8) 


that balances interlayer coupling energy against in¬ 
tralayer entropy. 

Additional interactions are expected to restrict the 
continuous symmetry to the discrete one {9 = 2Tmj3 with 
integer n) of the microscopic model. Of these locking 
terms, the most important are the RG-relevant ones that 

/oN 

link nearby layers, and the dominant one, cos Sz "^(r), 
acts between layers with separation three. All lock¬ 
ing terms that act between nearest-neighbour layers are 
RG-irrelevant; the leading example is KiKd^bzir))"^ — 
{dyhzir))"^] C0S(5i^^(r) -t 2{d,,hzir)){dyhz{r)) smSi^\r). 

We compute the RG flow m of the interlayer cou¬ 
plings K± and and the vortex fugacity y as a func¬ 
tion of i near the Gaussian fixed point of the uncoupled 
multilayer system. To leading order in the couplings 


d\nK^ 


ain£ 

97^3 


and 


dy 

ain£ 



1?>KJ 


K. 2- 




2 - 


i8a: 


( 9 ) 


Higher order contributions generate (initially zero in 
a model with only nn interactions) from K± and ATi, 
and also renormalise K, but are otherwise unimportant 
at small coupling. 

This description of the RG flow breaks down at the 
scale where the largest of Kj_, and ?/ is 0(1). If y ^ 1 
with K±,K 2 1, the RG treatment of weakly coupled 

layers can be applied at all scales: for y > 1, AT flows to 
zero, and hence so do K± and K^, yielding a conventional 
paramagnet. Alternatively, when K± ^ 1 with j/ <C 1, 
layers are strongly coupled at large scales and the state 
of the system depends on a competition between vortex 
and locking effects, which we now examine. 

To discuss the system of strongly coupled layers, we 
consider the energy cost of small amplitude fluctuations 
about a ground state of "Hsd- Let 


Then at quadratic order 


'Hsu = 


K 


2(27r)3 


d'q£(q)l<p(q)l' 


( 12 ) 


with (taking 9 = 0) 


£:(q) ^ql + (qy + K±smqz)'^ + «!(!- cosg^,)^ . (13) 

This dispersion relation is unusually soft (cx gf) in the 
interlayer direction. (The fluctuation energy for smectic 
liquid crystals [ 7 ] has a similar form, but with two quartic 
directions and a single quadratic one.) Because 
computed Eq. © is finite, the scaling flow of Eq. (§ 
stops at K± ~ 1. Without vortices, the system has long- 
range order and any non-zero locking interaction pins 9. 

To understand the influence of vortices on the system 
at K± > 1 we should examine ground states of Hsu hr 
height-held sectors with hxed vortex locations. The out¬ 
comes are, hrst, conhnement of vortex-antivortex pairs, 
and second, destruction of long-range order by bound 
pairs if Ks is small. 

Gonsider introducing a single vortex-antivortex pair in 


layer z = 0. Then, using the notation of Eq. (10), :/5o(r) 


is multiple-valued, winding by 6 (—6) along any closed 
path encircling the vortex (anti-vortex). At large vortex 
separation i^o(r) has a step of height 6 and width w on 
the line joining the vortex centres. The energy cost per 
unit length of this step is 


e ~ Kw{w ^ -I- K^) oc Kk_i 


(14) 


where the last expression follows from minimising over 
w and the optimal width is w ^ . Hence vortices at 

separations large compared to are subject to a linear 
conhning potential, and pairs are tightly bound. 

Tightly-bound pairs generate distortions in hz{x) that 
fall off only slowly with distance from the pair centre. 
The appropriate far held can be induced without consid¬ 
ering a multiply-valued v^o(r) by instead adding a cou¬ 
pling —K f d^r u(r)(/3o(r) to Rsu- We hnd [16] that the 
Fourier transform of the required potential w(r) has the 
form 


u(q) = 6i z • (q X b) (15) 


in the limit q k±, for a pair that has separation vector 
b and its centre at the origin. In the minimum-energy 
state containing this pair, the single-valued far held at 
r,z is (with q_L = {qx,qy)) 




1 

(27r)3 


d^c 


u(q±r+g,z) 

£(q) 


(16) 


hz{r) = hW(r) -kv5z(r) 
and write (with qj^ = (g^,, qy)) 

1 


( 10 ) 


Mr) = ji^J cl3q7>(q)e*^'»-'-+^^^) ■ (11) 


Extending this calculation to randomly located pairs 
at density p, we hnd that they generate huctuations in 
<Pz(j) with mean square amplitude 


{[MM) = 


(27r)3 


d^q 


k(q)P 

EM) 


(17) 
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The integral on the right of Eq. ( [T7| is divergent at small 
q, indicating that any non-zero density of bound vortex 
pairs destroys long-range order in the absence of lock¬ 
ing interactions. The in-plane and interlayer correlation 
lengths and can be estimated by using finite system 
size as a cut-off. Taking (|bp) ~ , we obtain the highly 

anisotropic results ~ and ^x ~ 

Spin correlations at wavevectors that have in-plane 
components close to one of the corners K of the trian¬ 
gular lattice Brillouin zone can be expressed in terms of 
small-qx components of hz{r). From Eq. ^ we find [TB] 

5(q) ^ (18) 

z 

where, on the right-hand side, the choice of ± sign and 
the value of the integer p depend on the zone corner. 

Correlations on scales shorter than ^z, '?_l resemble 
those in ground states and can be calculated from hi^\r) 
[Eq. ([^] by averaging over 6 . The result of this approxi¬ 
mation is a sharply defined helix 

TT Qiirp 7T 

S'(q) cx S{qj;-K±-cosQz) 6{qy-\ — ^±kx ^ sing^), (19) 

which is broadened when finite ^z, ^_l are taken into ac¬ 
count. 

Locking interactions suppress fluctuations in pz (r) and 
stabilise long range order if their effect integrated over 
the correlation volume is large at the RG scale on which 
K± ^ 1. The condition t)'^^z ~ 1 implies a crit¬ 

ical value of Kl ~ (p^^)®. Varying temperature in the 
spin system, a regime with strong interlayer correlations 
separates the high-temperature paramagnet from the 3D 
ordered phase if the renormalised K^, at the scale 

for which iLx ~ 1 - 

Since, from Eq. grows faster under RG than 

i^x, the existence of an intermediate regime requires 
both a sufficiently small microscopic value of J 3 and suf¬ 
ficiently slow generation of from K± and the RG- 
irrelevant coupling Ki. Because couples sites three 
layers apart, its value after an 0{1) RG rescaling in a 
system that, at the microscopic level, has only nearest- 
neighbour interactions cannot be larger than C>(Jx/T)” 
with n = 3. We find however (as noted in Ref. El for a 
similarly frustrated 2D problem) that this leading-order 
contribution is absent and the lowest non-vanishing con¬ 
tribution has n = 7. 

To obtain a phase diagram for the system, we in¬ 
tegrate the RG flow equations ([^ from a microscopic 
scale £0 to the scale i± at which K±{£±) = 1, with ini¬ 
tial values K± = J±/T, ~ {J^/TY and y{£o) = 

, The crossover between the conventional, high- 
temperature paramagnet and the cooperative paramag¬ 
net is at y{£i_) ^ 1, implying Jx ~ T[?/(f'o)]The 


ordering transition is at K 3 {£±_) = K^; setting = 
[vY)] in the expression derived above for this im¬ 
plies Jx T[y(£o)]^^^^. The cooperative paramagnet 
therefore occupies the range of interlayer couplings 

T[yYo)]^^^ <J±< T[y{£o)f^^ . ( 20 ) 

In the small J±/J limit of interest, this simplifies: the on¬ 
set of strong interlayer correlations is at Jj_ « 
while the ordering transition is at the much smaller value 
Jx ~ Gorrelations lengths in the coopera¬ 
tive paramagnet are given by ^z ~ Jx/J)"^ and 

C-L ~ ^o{J/Jj-Y^z- Hence the interlayer correlation 
length increases from « 1 at onset to ^z ~ 
near the ordering transition, and the in-layer correlation 
length from ^x ~ £o{J/J±Y to ^x ~ £oiJ/J±Y^^^- The 
parametrically large values reached by the correlation 
lengths are confirmation that the system indeed behaves 
as a cooperative paramagnet. 

In summary, we have discussed the nearest-neighbour 
Ising antiferromagnet on a lattice of weakly-coupled tri¬ 
angular layers with a frustrated abc stacking, showing 
that it has a low-temperature regime in which there 
are strong interlayer correlations without long-range or¬ 
der. This cooperative paramagnet is striking both for 
its correlations, which generate maxima in S'(q) on he¬ 
lices in reciprocal space, and for the mechanism of order- 
suppression via vortex-antivortex bound pairs, which is 
unlike that in any other system we are aware of. We 
expect that weakly coupled triangular layers in an hep 
lattice (the abab stacking) [THl [Hj will exhibit a similar 
cooperative paramagnetic regime. This system, together 
with a numerical study of the model with abc stacking, 
will be discussed in detail in forthcoming work m- 

An Ising model similar to Eq. has been proposed 
pni [ 2 T] as a description of charge ordering in the ma¬ 
terials LuFe 204 and YbFe 204 : Ising pseudospins rep¬ 
resent the charge states Fe^’*' and Fe^’*', and antiferro¬ 
magnetic coupling arises from Coulomb repulsion. Ex¬ 
perimental studies pniEUHn], in particular of YbFe 204 
[ 22 ) . find helices of scattering intensity in a temperature 
range above a three-dimensional charge-ordering transi¬ 
tion. While an accurate description of these materials 
would require treating both a more complicated (bi-layer) 
structure and additional (spin) degrees of freedom, the 
results we present in this paper elucidate how strong in¬ 
terlayer correlations can arise without long-range order. 
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dmr 1352271, JTC is supported in part by EPSRC Grant 
No. EP/I032487/1. 


[1] L. Balents, Nature 464, 199 (2010). 




5 


[2] G. Wannier, Phys. Rev. 79, 357 (1950); R. M. F. Houtap- 
pel, Physica 16 425 (1950). 

[3] S. T. Bramwell and M. J. P. Gingras, Science 294, 1495 

( 2001 ). 

[4] H. W. J. Blote and H. J. Hilhorst, J. Phys. A 15, L631 
(1982); B. Nienhuis, H. J. Hilhorst, and H. W. Blote, ibid. 
17, 3559 (1984); B. Nienhuis, Phys. Rev. Lett. 49, 1062 
(1982). 

[5] S. V. Isakov, K. Gregor, R. Moessner and S. L. Sondhi, 
Phys. Rev. Lett. 93, 167204 (2004) 

[6] M. K. Phani, J. L. Lebowitz, and M. H. Kalos, Phys. Rev. 
B 21, 4027 (1980); T. L. Polgreen, ibid. 29, 1468 (1984). 

[7] P. M. Chaiken and T. G. Lubensky, Principles of Con¬ 
densed Matter Physics (Gambridge University Press, 
1995). 

[8] O. A. Starykh and L. Balents, Phys. Rev. Lett. 98, 077205 
(2007). 

[9] D. Bergman, J. Alicea, E. Gull, S. Trebst and L. Balents, 
Nature Phys. 3, 487 (2007). 

[10] E. Rastelli, A. Reatto and A. Tassi, J. Phys C 16, L331 
(1983). 

[11] G. Zeng and C. L. Henley, Phys. Rev. B 55, 14 935 
(1997). 

[12] A. Danielian, Phys. Rev. Lett. 6, 670 (1961). 

[13] N. D. Mackenzie and A. P. Young, J. Phys. C 14, 3927 
(1981). 

[14] B. Nienhuis in Phase Transitions and Critical Phenom¬ 
ena edited by G. Domb and J. L. Lebowitz (Academic 
Press, London, 1987), Vol. 11. 

[15] J. Stephenson, J. Math. Phys. 11, 413 (1970). 

[16] See supplementary material. 

[17] D. T. Liu, F. J. Burnell, L. D. C. Jaubert and J. T. 
Ghalker, in preparation. 

[18] D. Auerbach, E. Domany, and J. E Gubernatis, Phys. 
Rev. B 37, 1719 (1988). 

[19] D. S. Zimmerman, G. Kallin, and A. J. Berlinsky, Phys. 
Rev. B 37, 7766 (1988). 

[20] Y. Yamada, S. Nohdo, and N. Ikeda, J. Phys. Soc. Japan, 
66, 3733 (1997); Y. Yamada, K. Kitsuda, S. Nohdo, and 
N. Ikeda, Phys. Rev. B62, 12167 (2000). 

[21] A. B. Harris and T. Yildirim, Phys. Rev. B 81, 134417 

( 2010 ). 

[22] A. J. Hearmon et al, Phys. Rev. B 85, 014115 (2012). 

[23] For a recent review, see: N. Ikeda, T. Nagata, J. Kano, 
and S. Mori, J. Phys Cond. Matt. 27, 053201 (2015). 


Supplementary material 

We provide derivations of the inter-layer coupling in 
terms of height fields, the far-held potential introduced 
by bound vortex pairs, and the structure factor in the 
cooperative paramagnet regime. 


HEIGHT MAPPING AND INTER-LAYER 
GOUPLING 

As hrst observed in Ref. [TJ the possible ground state 
conhgurations in a triangular lattice Ising antiferromag- 
net can be represented by a height held. The height 
representation we use here follows Ref. [H in which a 
height is assigned to the centre of each triangle. To do 
so, we choose a convention for dividing the triangular 
lattice into three sublattices A, B, and C. The six pos¬ 
sible ground-state spin conhgurations on a given triangle 
can then be identihed with 6 physically distinct height 
variables as follows: 


h (cta, ctb, CTc) 

0 (+> “> “) 

2 (-) -J +) 

4 (-, +, -) 


h ( cta , ctb , crc) 

1 (+) “> +) 

3 (—, +, +) 

5 (-b, -b, —) (21) 


The height held is single-valued provided that each tri¬ 
angle contains no more than one frustrated bond. Spin 
conhgurations containing triangles in which spins are all 
up or all down correspond to vortices in the height held. 

To obtain the inter-layer couplings, we work with 
the inverse mapping, given by ctq = fih -b Sq,), with 
(sa,sb,sc) = (0 ,2,-2), and 


fih) 


4 Tih 1 
3 cos — - - cos TTh 


( 22 ) 


Provided h is an integer, this reproduces exactly the spin 
conhgurations in Eq. (21), irrespective of which of the 6 
triangles neighbouring the site in question is used. In or¬ 
der to coarse-grain the system, we will use the convention 
CTz.a(i’) = Hbzir + a„) -b s„), where 


“a) = (s- = (^5' 

(23) 

On the a6c-stacked triangular lattice, the inter-layer 
Ising interaction couples each spin to three spins in each 
adjacent layer (one from each sublattice). In terms of the 
height helds, this gives: 


nA. = j± E E fihzir -b a^) -b Sa)/(r, z -b 1) 

a=A,B,C r 


( 24 ) 
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where 

7(i’,2)= 7(/» z+l (r + a„) + s„) (25) 

a=A,B,C 

Since we are concerned with the interlayer coupling only 
for the long-wavelength components of the height fields, 
we may Taylor expand this expression in derivatives of 
hz(r). The leading-order term is costt/ij, costt/i^+i which 
is highly irrelevant for weak inter-layer couplings. The 
only RG-relevant contribution, which comes from the 
first derivatives of h, is 

■H-L = (cos 7(h^+i(r) - h^{r))d^h^{r) 

^ r ,2 

TT \ 

-sin-(/iz+i(r) - h^{r))dyh^{r) j . 

(26) 

Approximating the sum by an integral (and omitting the 
factor 47r/9-\/3 for simplicity) reproduces the expression 
given in the main text. 


EFFECTIVE POTENTIAL DUE TO BOUND 
VORTEX PAIRS 


A key aspect of our argument for the existence of a 
cooperative paramagnetic regime is the fact that, due to 
the unusually soft dispersion along the stacking direction, 
even tightly bound vortex pairs are sufficient to frustrate 
long-ranged order in our model. Here we derive the form 
of the effective potential u(q) used to demonstrate this 
fact in the main text. 

We consider the far-held potential from a vortex- 
antivortex pair in an isolated layer. The vortex and an¬ 
tivortex are separated by a vector b. We take the pair to 
be centered at the origin. This pair is described by the 
height held conhguration 


h{x,y) 


3 

TT 


arctan 


/X + b/2 • 

\y + h/2-y) 


For |r| 3> |b|, we have: 


— arctan 


f X — b/2 • 
\y-h/2-y) 
(27) 


h{x,y) 


3 5 • (b X r) 

TT 


(28) 


or, equivalently, 

h q « 6z—- (29 

q 

The far-held height conhguration in Eq. ( |2^ can 
be induced by adding a potential term i'(q) to the 
effective Hamiltonian for the height held. Specih- 
cally, for an isolated layer, the effective Hamiltonian 
(A:/[27r]2)/d2q j^^s the 

minimum energy conhguration 


<7>(q)^ 


^z(q) 

^(q) 


^^z(q) 


Thus choosing a potential 

z;(q) = 6 iz • (q x b) (31) 

we recover the correct far-held conhguration. As shown 
in the main text, this potential term destroys the long- 
ranged order in ip. 


CORRELATION FUNCTIONS FROM HEIGHT 
VARIABLES 


Here we derive the expression for the spin-spin correla¬ 
tion function in terms of the height helds. Let the primi¬ 
tive lattice vectors within a triangular layer be = (1, 0) 
and a .2 = (1/2, ^/Z/2). Then the primitive lattice vectors 
of the afec-stacked triangular lattice have the form 

/2 1 

ai = (a'l, 0) , a2 = (a^, 0) , ag = ( -a '2 - -a^, 1 

(32) 

where we have taken the inter-layer distance to be 1 for 
convenience. Let bi, b 2 and ba denote the corresponding 
reciprocal lattice vectors. 

We hrst observe that the corners of the triangular lat¬ 
tice Brillouin zone lie at 

K = -bi -I- -b2 , K' = -bi -I- -b2 (33) 

Therefore, for a given position r = riiai + n 2 a 2 + on 
the stacked triangular lattice, we have 

47r Stt 

(K -H hihi -K /i 2 b 2 ) ■r=—{ni+ h 2 z) + —{n 2 - hiz) 

QiTt 47r 

(K' -H hihi + /i 2 b 2 ) • r = -^(ni - hiz) + —{n 2 + h 2 z) 

(34) 


modulo 27r. The combination ^(2ni -|- 712 ) mod 27r is 
simply Therefore 


gi(K+/iibi+/t2b2)'r _ gZ7r(So,-|-2pz)/3 
gi(K'-|-/iibi+/t2b2)'r _ —Sa+2pz)jZ 


(35) 


where the integer value of p = ( 2/12 — hi) depends on 
the Brillouin zone corner under consideration. Now, the 
most relevant terms in the spin-to-height mapping are 


^ ( 


gi 7 r /3 /i(r)gi 7 r/ 3 So 


g-iir/3 h(r)g-i7i-/3s, 


') • (36) 


(30) Keeping only slowly varying terms, we therefore obtain 
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(h(r)-h(0))^iz( 


^P+<lz 


g*qr 


(37) 


which gives Eq. (18) of the main text. 
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